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Abstract 

We  introduce  the  notion  of  a  discrete-time,  multivariate  point  process  which  can  arise  in 
the  modeling  of  an  optical  communication  system.  We  wish  to  estimate  the  rate  of  this  pro¬ 
cess  at  time  t  given  the  past  of  the  process  up  to  time  /- 1.  This  requires  the  computation  of 
a  certain  conditional  expectation:  we  perform  this  computation  by  introducing  an  absolutely 
continuous  change  of  measure  and  then  applying  the  generalized  Bayes'  rule.. 
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I.  Introduction 

Suppose  that  a  laser  beam,  whose  intensity  is  modulated  by  an  information  source, 
strikes  a  photodetector  which  is  a  part  of  the  receiver  of  an  optical  communication  system. 
Suppose  that  the  beam  is  also  influenced  by  external  random  phenomena;  this  occurs,  for 
example,  if  the  beam  passes  through  a  turbulent  atmosphere  [l].  The  occurrence  of  photoelec¬ 
trons  at  the  detector  can  then  be  modeled  as  a  doubly-stochastic,  time-space  Poisson  process 
[l].  It  is  often  necessary  to  compute  the  conditional  expectation  of  the  rate  of  such  a  process, 
for  example,  in  performing  a  likelihood  ratio  test  (Snyder  |2j,  Chapter  2  and  Chapter  6).  In 
general,  the  computation  of  this  conditional  expectation  is  quite  difficult,  although  results  are 
available  if  the  rate  process  has  a  Gaussian  form,  and  the  photodetector  surface  is  very  large; 
in  addition,  results  for  linear  estimates  are  available  when  the  photodetector  surface  is  arbi¬ 
trary,  but  the  rate  process  still  has  a  Gaussian  form  [3],  (see  also  [4,  5]  for  related  filtering 
results). 

Let  D  represent  the  photosensitive  surface  of  the  photodetector.  It  is  well  known  that 
the  probability  is  negligible  that  more  than  one  photoelectron  will  occur  in  the  entire  region 
D  during  a  time  interval  of  length  At  ,  if  At  is  sufficiently  small.  Let  DK  ,  be  a 

partition  of  D  into  disjoint  subregions.  Let  At  be  “sufficiently  small,”  and  let  n,  (k  )  denote 
the  number  of  photoelectrons  occurring  in  the  region  D *  ,  during  an  interval  (l ,  l  +  At  ].  We 
should  expect  that  for  each  1  <k  <I\  ,  n,(k)  takes  only  the  values  0  or  1.  In  addition,  we 
can  have  n,  (k  )  —  1  for  at  most  one  k  :  for  i  =£  k  .  we  must  have  n,  («  )  =  0.  In  this  paper  we 
formulate  the  discrete-time  version  of  the  rate-estimation  problem  by  introducing  the  notion 
of  a  discrete-time,  multivariate  point  process.  By  adapting  the  continuous-time  procedures 
found  in  (Bremaud  [6;,  Chapter  rvr;  see  also  pp.  69-70,  80  of  Chapter  III),  we  solve  the  prob¬ 
lem  of  estimating  the  rate  of  a  discrete-time,  multivariate  point  process.  (For  more  general 
discrete-time  procedures,  see  [7,  8]  ). 
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II.  Problem  Statement 

Let  K  >1  be  a  fixed  integer.  Let  (  fi,  F,  P  )  be  a  probability  space  on  which  all  the 
random  variables  in  this  section  are  defined.  We  call  {  n,  (A  );  t  >1,  l <k  <K  }  a  discrete¬ 
time,  multivariate  point  process  if  each  n(  (A- )  takes  only  the  values  0  and  1  and  if  the  events 
{n,  (A  )  =  1 } ,  1<A  <K ,  are  disjoint,  so  that  simultaneous  events  do  not  occur.  Next,  let  G0 
denote  the  trivial  cr-field  on  O,  and  set 

Gt  t  a{  n,(A);  l<s  <t ,  l<k  <K  }. 

Now,  let  X  be  a  sub-cr-field  of  F,  and  let 

Ft  t  Gt  v  X  ;  /  >0, 

denote  the  smallest  cr-field  containing  G  t  X .  We  assume  that  the  random  variables 

X,(A)  t  E  (  n,  (A )  j  F,_,  ]  ;  t  >1. 

are  actually  X-measurable  (note  that  X  =  F0)  with  known  joint  distributions.  For  example, 
we  might  take  X=<r(X')  for  some  random  variable  X  with  known  distribution  function 
F(x).  Then  each  X,  (A )  would  be  some  Borel  function  of  X .  Clearly,  the  joint  distributions 
of  the  {  X|  (A  );  t  >1,  1<A  </v  }  would  be  known,  at  least  in  principle.  Now,  observe  that 
X,  (A  )  =  P(  n(  (A  )  =  1  |  F,  _j  );  since  we  assume  that  simultaneous  events  do  not  occur, 

E  X.(A)  =  P(  U  {  n,(A)  =  i  }  |  F,_, 

k =1  k = I 

K 

This  implies  that  p,  =  1  -  E  Xt(A)  is  also  a  conditional  probability,  and  hence,  nonnega- 

* 

tive  and  bounded  above  by  1. 

Our  objective  is  to  compute 

X«(A)  *  E  [  n,  (A )  |  ] 

in  terms  of  the  known  joint  distribution  of  {  X,  (A  );  i<«  <f ,  1<A  <K  }. 
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III.  A  Reformulation  of  the  Problem 

We  shall  solve  our  problem  by  reformulating  the  probabilistic  setting  above  on  a 
different  probability  space,  (  n,  F ,  Q  ).  On  (  fi,  F ,  Q  )  let  {  ut  (k );  t  >1,  l<k  <K  }  be  a 
discrete-time,  multivariate  point  process.  In  a  manner  analogous  to  that  outlined  in  the  pre¬ 
vious  section,  we  take  C0  to  be  the  trivial  cr-field  on  fi,  and  set 

€,  ±  i<s  <t .  l<k  <K  }  ;  f  >1. 

Next,  let  X  be  a  sub-cT-fleld  of  F ,  and  set 

F ,  =  € t  v  X  :  l  >0. 


We  now  assume  that 

Eq  !  "i  (k  )  |  F ]  ==  pt,  {k  ),  (1) 

where  the  {/i,(A  )  }  are  arbitrary  constants  satisfying  pi,  (k )  >  o  and 
K 

qt  =  l  -  JD  /i(  (A- )  >  0-  (The  symbol  Eq  denotes  expectation  with  respect  to  the  measure 

k  =1 

Q).  As  a  consequence  of  the  above  assumption,  under  the  measure  Q,  each  i>,  ( k  )  is  indepen¬ 
dent  of  €,_!•  In  addition,  the  cr-fleld  G,  is  independent  of  the  <7-fleld  X  (see  Appendix). 
Next,  let  the  random  variables  {  X,  (k  ):  t  >  1,  1  <  k  <K  }  defined  on  fi  be  X-measurable 
and  have  the  same  joint  distributions  under  Q  as  {  t  >1,  l<)t  <K  }  (defined  on  fl  ) 

under  P. 

Given  the  preceding  probabilistic  setting  on  (fl,  F,Q),  we  make  the  following 
definitions.  Let 


vt  -  i-  E  M*)  :  1  >l 

k  —1 


A  K  (  \t{k)  p,  \ 

h  t  £  — 777  -  —  U(*>  + 

kZx  l  M*)  I 


p  i 

it 


t  >i. 


and 


(2) 
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Lt  t  n  i.  ■  <>i. 


with  L0  =  1.  Obsen'e  that  the  denominators  in  (2)  are  nonzero  and  deterministic.  Also, 
0<  X,  (k  )  <1  Q-a.s.,  and  hence,  p,  is  clearly  bounded  and  Q-integrable.  (Recall  that  the 


joint  distributions  of  the  {  X,  (k  )  }  under  Q  are  the  same  as  those  of  the  {  Xt  (Jt )  }  under  P). 

K  -  -  K 

Consequently,  l,  and  L,  are  Q-integrable.  Since  X,  (k )  =  l-pt,  and  y)  fit  (k ) 

k  =i  *  =  i 

=  1  -  q,  ,  it  is  easy  to  see  that  Eq  [  /,  |  i7,.,  ]  =  1.  Since  L,  =  ,  it  is  clear  that 

Eq  [  L,  |  IF  <_!  ]  =  Lt  ;  (3) 

i.e.,  L(  is  an  F  t  -martingale  under  Q.  Since.  Eq  [  L,  ]  —  Eq  [  L,  ]  =  Eq  [  /  ,  ]  =  1,  we  can 
define  a  new  measure  P  on  F ,  by 


P(F)  1  Jr  L,d Q  ;  Fe  F t  . 

(Technically,  we  should  show  that  the  family  {  Lt  }  is  uniformly  integrable.  However,  for 
our  purposes  this  is  not  necessary  since  if  we  wish  to  compute  E  [  |  £>_,  ]  for  some  r, 

we  can  select  any  finite  T  >t  and  then  restrict  our  attention  to  Observe  that 

since  L  0  =  1,  P  =  Q  on  JT .  If  E  denotes  expectation  with  respect  to  the  measure  P,  it  is 
not  hard  to  verify  (since  simultaneous  events  do  not  occur)  that 

E  {",(*)  I  E(.,]=EQ  [*/,(*)  |  Fh  i 

Lt- 1 

=  Eq  [  I 't(k)  I,  |  Fm]  (4) 

=  X,  {k  ), 

where  the  first  equality  in  (4)  follows  from  the  generalized  Bayes'  rule.  Now,  since  P  —  Q  on 

X, 

M*)  =  E[M*)  I  F,_,] 

under  P  is  probabilistically  equivalent  to  X,(A  )  under  P.  In  fact,  we  can  make  the  following 
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statement.  The  probabilistic  relationships  among 

{  nt  (A )  },  {  X,  (A )  },  {  Ft  },  {  Gt  },  and  X 

under  P  are  the  same  as  those  among 

{  vt  (A  )  },  {  X,  (A- )  },  {  F ,  },  {  G  t  },  and  X 

under  P.  In  the  next  section  we  shall  compute  E  [  vt  (A  )  |  G,_x  )  explicitly  as  a  Borel  func¬ 
tion  of  {  v,  (»' );  l<s  <t-l,  1<»  <K  }.  From  the  statement  above,  it  follows  that 

E  [n,(A)  |  ©(_,]  will  be  equal  to  the  same  Borel  function  applied  to 

{  n,  (t  );  l<«  <<-l,  i<t  <A"  }. 


IV.  Calculations 

In  this  section  we  compute  E  (  tst  (A  )  |  j.  By  Bayes’  rule, 

Eq  \ut(k  )Lt  |  Gt_x  ] 


E  [ut(k)  |  Gt_x 


EQ[Lt  | 


Next,  observe  that 


Eq  I  L(  |  G  I  _1  ]  —  Eq  (  Eq  [  Lt  \  F  t  ]  |  G  |  )  ■ —  Eq  I  Lt  |  G  ( . 

by  equation  (3).  In  addition, 

Eq  !  t'f  (A )  L,  |  G  ( ]  =  Eq  j  Eq  j  (A  )  L,  |  G t  ]  |  <!>(_,] 

—  Eq  j  Vt  (A  )  Eq  [  L,  |  G  t  }  |  C4_,  ]. 


Hence,  (5)  can  be  rewritten  as 


E  1  v,{k)  |  Gt_x 


Eq  1  M*)Eq  \Lt  |  G,  ]  t  Ct_x 


(5) 


(6) 


Eq  1  E|_i  j  G,_x  ) 

Before  proceeding  further,  it  will  be  convenient  to  introduce  the  following  notation.  For  each 
t  >1,  let 


t't  =  [  (1).  ■  •  •  ,  t'f  (K)  ]' 


where'  denotes  transpose.  Then  ut  is  a  K  -dimensional  random  vector.  Now,  let 


-  (  "i.  •  -  •  ,  "t  1'  • 


Clearly,  ut  is  a  l/\T-dimensional  random  vector.  Next,  let 

*«  =  [  *(»),  ■  •  ■  .  -,(K)  ]'  and  7,  t  [  _ zt  ]' 

denote  “dummy”  variables  in  RA  and  JR(A  ,  respectively.  Define 


A'  X,  (A )  pt  p. 

M*)  =  £  —777  -  —  *.(*)  +  — , 

/'«(*)  9t  <?< 


(7) 


and 


htUt)  = 


[  X,(t) 
I  h,  (A  ) 


V  t 
Qt 


,,  .  V  t 

2t(k)  +  -  • 

Qt 


Clearly,  ht(zt )  is  an  JF-measurable  random  variable  and  h  ,(z, )  is  an  AT-measurable  random 
variable.  In  fact,  h  t(zt)  under  Q  is  probabilistically  equivalent  to  ht  (zt  )  under  P.  It  is  then 
clear  that  the  same  can  be  said  of 


H,(zt)  t  H  h,(z,), 
•  =1 

and 


(8) 


Ht(jt)  t  n  mo. 

$  =1 

Next  observe  that  lt  —  A  ,  (iy,  ),  and  that  L,  =  We  are  now  ready  to  compute 

Eq  (  Lt  |  €t  ]  in  (6)  for  each  t  >1.  (Note  that  Eq  [  L0  |  C0  ]  =  Eq  [  L0  ]  =  1). 
Observe  that  L(  =  H  t(Pt)  is  an  X-measurable  function  of  a  G ,  -measurable  random  vector. 
Since  A  and  C,  are  independent  under  Q,  it  follows  that 

Eq  [Ht(Vt)  I  C,  ]  =  EQ[tft(7,)]|-_- . 


From  the  remarks  preceding  equation  (8),  we  have 
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EQ  !L(  |  C,  ]  =  E  [  H,(It)  )  |  _  =  .  . 

We  now  set 

/,  (7, )  =  E  (  H,  (7,  )  j.  (9) 

Note  that  /,  (7, )  is  a  deterministic  function  of  7,  .  Equation  (6)  becomes 

E  (  i/(  (k  )  |  ]  =  Eq  1  1  C‘-'  1 

/  ( -lC^t  -i) 

We  can  write 

Eq  !  t>,  (k)ft  (P( )  |  <&,_,]=  Eq  [  u,(k)  / 1  (77,  .  ut)  | 

Now,  equation  (1)  implies  that  //,  is  independent  of  C ,_t  under  Q.  Therefore. 

Eq  [»/,(*)/,  (17,)  I  C*-,  )  =  Eg  [!/,(*)/,  I/,)]|T  p  - 

t  -1  v\  -1 

Since  simultaneous  events  do  not  occur,  it  is  trivial  to  compute 

Eq  [  t'.(k)  f  t  (7,  ,  i/,))  =  ft(Tt_x,  e*)Q(  u,  (k  )  =  l  ) 

=  (k  )  ft  (7,.„  ek  ), 

where  ek  is  the  standard  unit  vector  in  Kk  with  a  1  in  the  k  th  position  and  a  0  in  the  other 
K- 1  positions.  We  conclude  that 


E  [  ut(k)  | 


i,  ek  ) 
ft-iPt-i) 


and  hence. 


E  (  n,  (k  )  | 


ek) 

f  t  -i(nt  -i) 


with  the  obvious  meaning  of  the  symbol 
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V.  Summary 

We  have  shown  that  if  {  nt  (k  )\  t  >l,  l<k  <K  }  is  a  discrete-time,  multivariate  point 
process  residing  in  the  probabilistic  setting  outlined  in  Section  n.  then 


E  [  n,  (k  )  |  £',_i 


(k  )ft  (nt-i,  ek) 

/  t  _,(«£  _l) 


where  f  t{zt)  is  given  by  equations  (7),  (8),  and  (9). 

The  introduced  in  Section  III  were  arbitrary;  if  we  set  /i,(fc)  —  — - — 

A  +1 

and  all  k  in  the  preceding  equations,  we  find  that 


for  all  t 


E  i  n,  (k  )  | 


Ck  ) 

a,  _i(nt  _,) 


where 


at  (7t )  1  E  [  n 


E  (E  (*)  -  v,  K  (*)  +  p. 


i  U=i 


Appendix 

For  1<«  <f  ,  let  z ,  denote  either  the  zero  vector  in  IR*  or  any  standard  unit  vector  in 
.  In  this  appendix  we  prove  that  if  E  €  X ,  then  (using  the  notation  of  Sections  III  and 

IV) 

Q(  I't  =*<,•••,  vx  =  z E  )  =  |  n  Q(  )  |  Q(E  )• 

Since  C  t  —  <r{  i/l,  .  .  .  ,  u(  },  this  will  prove  that  C  (  and  X  are  independent  under  Q. 

Proof.  Using  the  definition  of  conditional  probability, 

Q(  I /,=  2T, .....  I/,  =  z„  E  )  =  /  Q(  =  ■«(  I  F«_,  )  <fQ.  (10) 

{  - “  *  |.  £  } 
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Since  simultaneous  events  do  not  occur,  if  zt  =  ek  .  Q(  ut  =  zt  \  F )  ==  /i,  (it ).  If  z,  =  0, 

K 

Q(  vt  —  zt  |  F  t )  —  1  -  y]  nt  (k  ).  Since  Q(  i/t  —  zt  |  JF,_,  )  is  deterministic, 

*  =i 

Q(  Vt  =  z,  |  F )  =  Q(  ut  =  zt  ). 

Hence  (10)  becomes 

Q(  "i  =  *t .  •  •  •  .  ui=  tx.E)  =  Q(  v,  ==  zt  )Q(  2,_,,  .  .  .  ,  «/,=  r,,  E  ). 

The  remainder  of  the  proof  by  induction  is  clear. 

QED 
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